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The  U.S.  Department  of  Energy's  (DOE)  Hydrogen  Storage  Engineering  Center  of  Excellence  (HSECoE)  is 
developing  models  for  hydrogen  storage  systems  for  fuel  cell-based  light  duty  vehicle  applications  for  a 
variety  of  promising  materials.  These  transient  models  simulate  the  performance  of  the  storage  system 
for  comparison  to  the  DOE's  Technical  Targets  and  a  set  of  four  drive  cycles.  PNNL  developed  models  to 
simulate  the  performance  and  suitability  of  slurry-based  chemical  hydrogen  storage  materials.  The 
storage  systems  of  both  a  representative  exothermic  system  based  on  ammonia  borane  and  an  endo¬ 
thermic  system  based  on  alane  were  developed  and  modeled  in  Simulink®.  Once  complete,  the  reactor 
and  radiator  components  of  the  model  were  validated  with  experimental  data.  The  system  design  pa¬ 
rameters  were  adjusted  to  allow  the  model  to  successfully  meet  a  highway  cycle,  an  aggressive  cycle,  a 
cold-start  cycle,  and  a  hot  drive  cycle.  Finally,  a  sensitivity  analysis  was  performed  to  identify  the  range 
of  material  properties  where  these  DOE  targets  and  drive  cycles  could  be  met.  Materials  with  a  heat  of 
reaction  >11  kj  mol-1  H2  generated  and  a  slurry  hydrogen  capacity  of  >11.4%  will  meet  the  on-board 
efficiency  and  gravimetric  capacity  targets,  respectively. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Engineering  of  onboard  hydrogen  storage  materials  for  fuel  cell 
applications  is  one  of  the  key  objectives  of  the  U.S.  Department  of 
Energy's  (DOE)  Hydrogen  Storage  Engineering  Center  of  Excellence 
(HSECoE)  program.  The  goal  is  to  model  and  optimize  the  necessary 
hardware  required  to  build  hydrogen  storage  systems  for  light  duty 
vehicle  applications  1  .  On  the  basis  of  their  promising  H2  storage 
capacity,  ammonia  borane  (AB)  and  alane  have  been  chosen  by 
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HSECoE  for  the  initial  chemical  hydride  simulation  studies.  Several 
studies  on  the  science  and  material  properties  of  AB  and  alane  have 
been  reported  in  the  literature.  Autrey  and  co-workers  studied 
hydrogen  release  from  AB,  its  thermal  stability,  and  other  engi¬ 
neering  aspects  [2,3],  while  Aardahl  et  al.  discuss  system  issues, 
such  as  heat  management  and  reactor  design,  that  need  to  be 
considered  for  vehicular  storage  system  designs  [4].  Several  re¬ 
searchers  have  studied  the  decomposition  kinetics  of  alane  and 
developed  kinetic  expressions  for  hydrogen  generation  [5-7]. 
Alane  also  has  been  studied  for  vehicle  applications  [8-10], 
although  a  significant  concern  is  its  regeneration  [11  .  Our  work  at 
HSECoE  focuses  on  designing,  developing,  and  evaluating  reactor 
systems.  We  use  a  modeling  approach  based  on  scientific  knowl¬ 
edge  for  chemical  hydrogen  storage  materials  reported  in  the 
literature. 


K.P.  Brooks  et  al.  /  Journal  of  Power  Sources  271  (2014)  504—515 


505 


While  there  have  been  a  large  number  of  models  for  metal  hy¬ 
dride  storage  systems  reported  in  the  literature,  Pacific  Northwest 
National  Laboratory  (PNNL)  was  the  first  to  develop  systems 
models  of  various  chemical  hydrogen  storage  materials  (solid  AB, 
solvated  AB,  and  alane)  for  engineering  analysis  and  development 
[12-16].  Chemical  hydrogen  storage  materials  differ  from  metal- 
hydride  storage  systems  in  that  they  must  be  removed  from  the 
vehicle  before  they  can  be  regenerated  by  chemical  processing. 
Metal-hydride  storage  materials  allow  direct  addition  of  hydrogen 
onboard  a  vehicle  for  regeneration. 

Previous  studies  using  PNNL-developed  systems  assumed  that 
the  chemical  hydrogen  storage  materials  were  a  pelletized  solid 
that  was  reacted  either  in  place  or  by  moving  through  the  system 
using  auger  or  pneumatic  transport  [12-14  .  Because  of  the 
challenge  of  moving  solid  pellets  through  a  pressurized  system 
while  preventing  hydrogen  loss  to  the  environment,  the  use  of 
solids  has  been  rejected  in  the  HSECoE.  Alternatively,  pure  liquid 
systems  would  be  ideal  in  terms  of  transport  both  within  and 
outside  the  vehicle.  Studies  also  have  been  performed  assuming 
AB  is  dissolved  in  an  ionic  liquid  [15,16].  However,  successfully 
developing  a  solvated  fluid  at  high  concentration  that  remains  a 
liquid  both  pre-  and  post-reaction  has  been  elusive.  As  a  result,  the 
focus  of  the  research  for  the  HSECoE  has  been  the  development  of 
these  materials  in  a  slurry  form.  This  ongoing  development  work 
has  demonstrated  AB  slurries  in  silicone  oil  at  solids  concentra¬ 
tions  up  to  45  wt%  [17  and  alane  in  either  silicone  oil  or  me¬ 
chanical  pump  fluid  at  solids  concentration  up  to  60  wt%.  This 
paper  focuses  on  the  development  of  a  system  and  a  model 
describing  the  transient  behavior  of  a  slurry-based  hydrogen 
storage  system.  The  model  will  include  both  AB  and  alane. 
Experimental  validation  of  these  models  has  been  performed  for 
two  key  system  components.  The  experimental  work  associated 
with  this  testing  will  be  described. 


2.  Slurry  system  design 

The  slurry  system  design  on  which  the  modeling  and  experi¬ 
mental  work  is  based  is  shown  in  Fig.  1.  Details  of  this  system  and 
how  components  were  selected  are  provided  elsewhere  [18],  but  a 
brief  description  of  the  process  will  be  provided  here.  In  this  sys¬ 
tem,  the  slurry  is  transported  on  and  off  the  vehicle  between  the 
filling  station  and  a  volume  displacement  tank.  A  volume 
displacement  tank  configuration  minimizes  the  system  compo¬ 
nents  by  requiring  only  one  tank  for  both  the  feed  and  the  product. 
A  membrane  separates  the  feed  side  of  the  tank  from  the  product. 
As  fresh  slurry  is  loaded  into  the  tank,  the  membrane  forces  spent 
slurry  out  of  the  tank  and  back  to  the  filling  station  for  eventual 
reprocessing. 

During  startup,  electrical  resistance  heaters  bring  the  reactor 
walls  to  reaction  temperature.  Meanwhile,  the  fresh  feed  in  the 
volume  displacement  tank  is  mixed  to  ensure  homogeneity  before 
it  is  sent  into  the  system.  Once  the  reactor  reaches  the  required 
temperature,  the  feed  pump  begins  moving  fresh  feed  into  the 
reactor.  As  the  slurry  flows  through  the  heated  tubular  reactor, 
hydrogen  is  generated  in  a  non-catalytic  thermolysis  process.  The 
hydrogen  gas  is  separated  from  the  three-phase  mixture  in  a  phase 
separator,  cooled  in  a  radiator,  and  then  sent  back  to  the  volume 
exchange  tank. 

For  the  exothermic  AB  slurry,  a  fraction  of  the  hot  slurry  product 
from  the  phase  separator  and/or  cold  product  after  the  radiator  can 
be  recycled  to  the  feed  if  needed.  These  two  recycle  streams  preheat 
the  feed  stream  to  reaction  temperatures  while  adding  thermal 
mass  to  minimize  the  temperature  excursions  associated  with  the 
exothermic  reaction.  For  the  endothermic  alane  slurry,  the  feed  is 


not  recycled.  Instead  it  passes  through  a  recuperator  to  recover  the 
heat  of  reaction  before  being  cooled  in  the  radiator. 

Once  the  gases  generated  in  the  reactor  are  separated  from  the 
slurry,  they  pass  through  a  coalescing  filter,  a  radiator,  and  an 
impurities-adsorbent  bed  to  condition  the  hydrogen  feed  prior  to 
the  fuel  cell.  The  hydrogen  must  be  cooled  with  a  radiator  to 
below  85  °C,  and  impurities  must  be  removed  to  the  levels 
required  by  SAE  J2719  and  ISO/PDTS  14687-2,  which  is  >99.97% 
purity  19].  Rather  than  heat  the  system  quickly  to  meet  the 
stringent  startup  time  required  by  the  DOE  technical  targets, 
hydrogen  gas  is  stored  in  a  ballast  tank.  The  ballast  tank  allows  for 
immediate  flow  of  hydrogen  to  the  fuel  cell  during  cold  startup 
while  the  storage  system  is  warming  to  the  required  temperature 
to  initiate  the  thermolysis  reaction.  The  ballast  tank  also  provides 
hydrogen  to  accommodate  the  transients  during  operation. 

3.  Slurry  system  model 

The  mass  and  energy  balances  of  the  entire  storage  system 
described  in  the  previous  section  were  modeled  using  the  Simu- 
link®  simulation  software,  which  provides  an  understanding  of  the 
operational  dynamics  of  the  storage  system  under  simulated  drive- 
cycle  scenarios.  The  components  that  are  modeled  are  described 
below. 


3.1.  Reactor  model 


The  reactor  for  both  the  AB  and  alane  slurry  is  a  tubular  reactor 
with  an  internal  stirrer  or  auger  to  prevent  agglomeration  of  the 
solids  during  reaction  while  allowing  the  gas  to  flow  through  the 
reactor  center.  The  parameters  used  for  the  reactor  model  are 
provided  in  Table  1.  For  the  purposes  of  modeling,  it  is  assumed  that 
the  slurry  is  pushed  to  the  walls  of  the  reactor  while  the  gas  flows 
down  the  reactor  center. 

The  fractional  conversion  of  the  chemical  hydrogen  storage 
material  (CH)  to  H2  in  the  slurry  through  the  length  of  the  reactor  is 
given  by  a  and  is  calculated  as  shown  in  Eq.  (1). 


dm  dm  dm 

gF  =  -Uslurryg^  +  -gF 


kinetics 


In  the  case  of  AB,  Eq.  (1 )  consists  of  two  equations  (1  =  2)  for  the 
two  separate  polymerization  reactions  that  correspond  approxi¬ 
mately  to  the  first  equivalent,  and  the  second  and  a  fraction  of  the 
third  equivalent  for  a  total  H2  production  of  2.35  equivalents.  In  the 
case  of  alane,  Eq.  (1)  is  a  single  equation  producing  1.5  equivalents. 
The  velocity  of  the  slurry  (usiUrry)  is  calculated  based  on  equation 
Eq.  (2): 


_  ^slurry  7 

“slurry  —  “ - 4 - r  ^auger^auger 

/?slurry/1slurry 


where  Zauger  is  the  axial  displacement  of  the  auger  per  rotation, 
Gauger  is  the  rotational  speed,  and  AsiUrry  is  the  area  of  the  gas/slurry 
interface  in  the  reactor  where  the  slurry  channel  is  estimated  to  be 
75%  of  the  total  area  of  the  reactor  calculated  from 
n(dfn-d^ urry)/4.  As  hydrogen  is  generated,  the  slurry  density, 
Psiurry*  decreases  to  create  a  foamy  mixture  of  solid,  liquid,  and  gas 
with  a  density  given  in  equation  Eq.  (3). 


_  “^slurry  3“ 

y  171^2  .  nislurry/inert  _i_  ^tlsluriyTcH  ( 1  ai)  .  ^slurry/cH  OL\ 

Ph2  Pinert  PCH  PCHP 

The  parameters  f[neYt  and  fc h  are  the  fractions  of  inert  and 
chemical  hydride;  and  ph2,  pinert,  Pc h»  and  Pchp  are  the  densities  of 
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Fig.  1.  Slurry  AB  system  process  flow  diagram. 
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hydrogen,  slurrying  agent,  chemical  hydride,  and  its  product, 
respectively.  A  minimum  value  of  the  slurry  density  (psiurry,min)  is 
set,  and  as  more  hydrogen  is  generated,  it  then  goes  into  the  gas 
phase  rather  than  the  slurry.  This  foaming  of  the  slurry  mixture  has 
been  observed  during  experiments. 

The  rate  of  change  of  a  due  to  kinetics  is  given  by  Eq.  (4), 
following  the  Kolmagorov-Johnson-Mehl-Avrami  (KJMA) 
approach  based  on  sigmoidal  conversions  as  a  function  of  time, 
characteristic  of  nucleation,  and  growth  phenomena  in  solids 

[4]. 


rH; 


^  slurry/3  slurry  ^ChAZ  Q  dcq 

2^  P*  aF 


MW, 


kinetics 


The  exothermic  heat  generated  or  required  from  the  reactions 
is  critical  to  the  calculation  of  reactor  temperature  and  is  calcu¬ 
lated  based  on  the  reaction  kinetics  as  shown  in  Eq.  (7). 


^slurry^CH  Pslurry 


MWch 


A  ^ 

kinetics 


=  njIfiCl  -  «i)[  -  ln(l  -  1)/n,  (4) 

kinetics 

ki  is  the  rate  constant  of  the  dehydrogenation  reactions  calculated 
by  the  Arrhenius  expression  as  shown  in  Eq.  (5). 

ki=A'ie~E°-,/RT  (5) 


9 off 

at 


An  energy  balance  across  the  reactor  is  used  to  calculate  the 
temperatures  of  the  slurry  in  the  reactor  (Treac tor),  reactor  wall 
(Tmetai).  and  hydrogen  gas  (TH2)  as  shown  in  Eq.  (8).  The  model 
assumes  convection  of  mass  and  heat  but  no  axial  conduction  of 
heat  or  diffusion  of  mass.  The  source  term  (Qioss)  is  the  total  heat 
loss  calculated  from  the  metal  to  the  ambient  through  the  reactor's 
insulation.  During  reactor  startup,  heaters  on  the  outer  surface  of 
the  reactor  walls  conduct  heat  into  the  reactor  (Qheater)- 


a  „  r  ( ^reactor 
^slurry Pslurry'-pslurry  l  ' 


d^reactoiA 

Uslurry - ^ - J  +  dout^ftmetal (^reactor  —  Tmetal)  +  ^in7r^H2  (^reactor  —  ^H2) 


Or 


xn 


AnetalPmetal Cp, metal  ^  ^  ^  F  dout^^meta\(Tmeta\  ^reactor)  —  Qheater  Qioss 

/97lj  9TYj  \ 

AH2PH2C^2  \  +  UgaS_a z)  +  din7zhH2  (TH2  -  ^reactor)  =  0 


(8) 


The  overall  hydrogen  production  rate  (rHz)  is  calculated  using 
the  expression  shown  in  Eq.  (6).  The  hydrogen  coefficients  (ft) 
correspond  to  the  molar  ratio  of  hydrogen  generated  in  each 
dehydrogenation  reaction,  xch  is  the  fraction  of  CH  in  the  slurry,  A z 
is  the  length  of  a  differential  reactor  section,  and  MWch  is  the 
molecular  weight  of  CH. 


Table  1 

Reactor  parameters  used  in  the  reactor  model. 


Value 

Units 

Reactor  parameters 

Reactor  outer  diameter  (dout) 

4.44  x  1(T2 

m 

Reactor  inner  diameter  (djn) 

4.11  x  10  2 

m 

Reactor  stirrer  diameter  (dstir) 

4.76  x  10~3 

m 

CH  fraction  in  slurry  (xCh) 

0.5 

wt  frac 

Start-up  heater  (Qheater) 

8 

kW  ITT1 

Density  metal  (pmetai) 

8000 

kg  ITT3 

Heat  capacity  metal  (Cp>metai) 

480 

J  kg-'  I< 

Slurry  thermal  conductivity  (k'CH) 

2.0 

W  m  1  K 

Silicone  oil  density 

1070 

kg  ITT3 

Silicone  oil  heat  capacity 

1846 

J  kg-'  I< 

AB  parameters 

AB  bulk  density  (pbuik) 

780 

kg  ITT3 

Molecular  weight  AB  (MWCh) 

30.8 

g  mol-1 

Heat  of  reaction  1st  Eq.  (A HO 

15 

kj  mol-1  H2 

Heat  of  reaction  2nd  Eq.  (A H2) 

22 

kj  moD1  H2 

Heat  capacity  slurry  (Cpt)uik) 

2694 

J  kg-'  I< 

Reactor  length  (AB) 

0.5 

m 

Alane  parameters 

Slurry  bulk  density  (pbuik) 

1486 

kg  ITT3 

Molecular  weight  alane  (MWCh) 

30 

g  mol-1 

Heat  of  reaction  (AHi) 

-7.6 

kj  mol-1  H2 

Heat  capacity  slurry  (Cpbuik) 

1000 

J  kg-'  K 

Reactor  length  (alane) 

2.0 

m 

CH  =  chemical  hydrogen. 


The  cross  sectional  areas  of  the  metal  and  gas  are  calculated 
based  their  corresponding  diameters.  The  heat  transfer  coefficients 
between  the  metal,  the  slurry,  and  the  gas  (fimetai/H2 /slurry)  are 
calculated  based  on  the  appropriate  thermal  conductivities,  di¬ 
ameters,  and  the  Nusselt  numbers.  The  Nusselt  number  for  the 
slurry  at  the  gas  and  metal  interface  were  fitted  from  graphs  pro¬ 
vided  for  fully  developed  laminar  flow  in  concentric  annular  ducts 
for  uniform  temperatures  at  both  surfaces  from  Rohsenow  et  al. 
[20]  as  shown  in  Equations  (Eq.  (9)— (11)).  The  Nusselt  number  for 
the  gas  phase  was  equal  to  4.36  when  laminar  and  Equation  (Eq. 


(11))  when  turbulent  [21  . 

Nuin  =  3.75r*-°'436 

(9) 

Nu0Ut  =  3  +  r* 

(10) 

NuHi  =  0.023Re/.8Pr/.4 

(ID 

4  mH, 

UU2=  ,  2  X 

(12) 

PH2  7T  ( 1 ^siurry  —  ^stir  ) 

The  value  of  r*  is  the  ratio  of  the  inner  and  outer  diameter  of  the 
slurry.  ReHl  and  PrH2  are  the  Reynolds  and  Prandtl  numbers, 
respectively,  for  the  hydrogen  gas  based  on  a  gas  velocity  calculated 
from  its  mass  flow  and  the  area  between  the  gas/slurry  interface 
and  the  stirrer. 

A  variable  step  solver  “ode23tb”,  available  in  Simulink,  is  used 
for  all  the  simulations  conducted  in  this  work.  To  avoid  algebraic 
errors  in  the  model  resulting  from  the  reactor  outputs  feeding  back 
to  the  controller,  unit  delays  are  used.  Simulink  solves  the  temporal 
components  of  the  equations.  The  spatial  components  are  solved  by 
a  series  of  10  differential  reactor  length  elements  coupled  together. 
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Time  (minutes) 

Fig.  2.  Batch  kinetic  results  for  a  35  wt%  AB  slurry  with  model  fits. 


3.2.  Kinetic  parameter  development 

The  Arrhenius  values  for  AB  were  fit  based  on  batch  reactor  data. 
Batch  reactor  data  of  35  wt%  AB  slurries  in  AR-20  silicone  were 
performed  by  the  researchers  at  the  University  of  Oregon  [22]. 
Hydrogen  production  as  a  function  of  time  over  the  range  of 
120-150  °C  is  shown  in  Fig.  2.  The  data  were  fit  to  obtain  the  rate 
constants  as  described  in  equation  Eq.  (5)  and  hydrogen  coefficients 


radiator  and  from  the  radiator  to  ambient.  The  radiator  is  assumed 
to  have  finned  tubes  with  ambient  air  flowing  across  them.  The 
temperature  distribution  in  the  fin  is  assumed  to  be  constant,  while 
the  base  metal  temperature  varies.  The  model  assumes  transient 
thermal  mass  of  the  radiator  with  convection  between  the  process 
stream  (rsiurry)  and  the  radiator  metal  (Tmetai)-  Axial  conduction  is 
considered  negligible.  The  equations  used  to  model  the  process  and 
metal  temperatures  are  shown  below  in  Eq.  (14). 


Anterior  Pslurry  ^-pslurry 


'61, 


slurry 


+  ^slurry  ^  +  Anterior^interior71  (^metal  ^slurry)  —  0 


dt 


AxteriorPmetal Cp, metal  f  ^  j  +  Anterior^interior'Tr (^metal  ^slurry)  +  Axterior^exterior^fin^fin  (Tmetai  ^ambient)  —  0 


(14) 


(ft)  described  in  equation  Eq.  (6)  that  would  be  used  in  the  Simulink 
reactor  model.  The  best  fit  with  this  experimental  data  was  ach¬ 
ieved  when  (3  was  found  to  be  a  function  of  temperature  rather  than 
constant  (Eq.  (13)): 


where  W  and  n!  are  constants  to  best  fit  the  hydrogen  coefficient's 
temperature  dependence.  The  Arrhenius  values  for  alane  were 
available  from  the  literature  and  were  based  on  those  previously 
developed  by  Graetz  and  Reilly  for  a-AlH3  polymorphs  [7  .  In  the 
case  of  alane,  the  (3  value  remained  constant  with  temperature.  As  a 
result,  no  exponential  fitting  was  required.  The  kinetic  and 
hydrogen  coefficient  parameters  for  the  AB  and  alane  slurry  model 
are  shown  in  Table  2. 

3.3.  Gas  and  slurry  radiators 

The  radiator  models  used  for  both  gas  and  slurry  product 
cooling  are  based  transient  heat  losses  from  the  process  to  the 


The  increased  surface  area  associated  with  the  fin  <5fin  can  be 


Table  2 

Reaction  parameters  developed  from  experimental  data. 


Parameter 

Units 

AB 

1st 

equivalent 

2nd/3rd 

equivalent 

Alane 

Arrhenius 
factor  (Ai) 

s”1 

9.6  x  103 

8.8  x  107 

1.2  x  1010 

Activation 
energy  (Ei) 

kj  mol-1 

45.0 

89.1 

102.2 

H2  coefficient 
constant  (/?•) 

mol  H2/mol  AB 

1.45 

0.905 

1.5 

Beta  constant  (/<• ) 

(1  /°C)ni' 

1.29  x  10  6 

1.62  x  10~7 

0 

Beta  temperature 
exponent  (n'f) 

Unitless 

2.98 

3.84 

1 

H2  coefficient 
exponent  (n,) 

Unitless 

3.1 

1.2 

2 
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calculated  by  Eq.  (15)  as: 


<5fi 


in 


nz  _  nz 
fin  exterior 


+  2DexteriorPfm  j  (^2DexteriorPfin^ 


(15) 


The  reactor  inlet  and  outlet  temperatures  are  controlled  by 
adjusting  the  fraction  of  spent  product  that  is  recycled. 

The  feed  flow  rate  is  based  on  the  hydrogen  requested  and  the 
hydrogen  capacity  of  the  storage  material  based  on  complete 
conversion  as  shown  in  Eq.  (18): 


where  Dfin  is  the  outer  diameter  of  the  fin,  Dexterior  is  the  outer 
diameter  of  the  heat  exchanger  tube,  TWoor  is  the  inner  diameter 
of  the  heat  exchanger  tube,  and  pfin  is  the  fin  pitch.  The  fin  effi¬ 
ciency  (rjfin)  was  calculated  for  a  circular  fin  as  described  by  Roh- 
senow  et  al.  [20].  The  heat  transfer  coefficient  for  a  cylinder  in  cross 
flow  is  calculated  using  on  a  method  reported  by  Incropera  and 
DeWitt  [21  as  shown  below  in  Eq.  (16): 


h 


exterior 


Kir 


0  62Re1/2  •  Prl‘. 

q  3  _j_  ^Dexterior- 1  'air 


1/3 


D 


exterior 


1  +  (0.4 /Prair)2/3 


1/4 


X 


1  + 


DP  \  5/8"]  4/5' 

^^Dexterior 


(16) 


ior  \ 


28200  ) 


where  Red  —  Pair^air^hxOD/Mair  ^nd  Pra ir,  /<air>  Hair*  Main  ^nd  Pair  ^1*6  the 
Prandtl  number,  thermal  conductivity,  velocity,  viscosity,  and  den¬ 
sity  of  air  at  the  film  temperature,  respectively.  The  value  of  ^interior  is 
calculated  based  on  the  Nusselt  number  for  the  interior  of  a  tube 
with  a  swirl  flow  device  for  laminar  flow  in  a  uniformly  heated  tube 
as  described  by  Rohsenow  et  al.  20]  and  shown  in  Eq.  (17): 

^interior  =  5.172  (l  +  0.005484*Pr°  7(ReDinterior/yf25) 1/2 

(17) 

where  y  is  the  twist  ratio  of  the  turbulator  inside  of  the  heat 


(18) 


where  nH2  required  is  the  molar  flow  of  hydrogen  requested  from  the 
fuel  cell.  To  maintain  the  system  pressure,  the  total  hydrogen 
required  is  adjusted  up  or  down  by  a  proportional-integral  (PI) 
controller.  The  ballast  tank  set  point  is  25  bar. 

The  flow  of  hot  slurry  product  as  recycle  with  the  feed  was 
adjusted  to  ensure  that  target  temperature  ( Ttarget)  was  met  and  the 
reaction  could  be  initiated  without  requiring  external  heating.  The 
flow  rate  of  this  recycle  stream  was  calculated  using  Eq.  (19). 


^recycle  — 


^feed^Pfeed  (^target  ^ambient) 
Q^product  (^product  —  ^target ^ 


(19) 


With  the  large  exothermic  heat  release  of  AB,  the  temperature  of 
the  product  exiting  the  reactor  in  some  cases  was  unacceptably 
high.  To  address  this  problem,  a  fraction  of  the  cold  product  from 
the  radiator  was  mixed  with  a  fraction  of  hot  product  exiting  the 
phase  separator  to  produce  an  overall  larger  recycle  stream  at  the 
same  target  temperature  as  before.  The  amount  of  recycle  flow 
from  the  radiator  and  phase  separator  is  calculated  using  Eq.  (20) 
and  Eq.  (21),  respectively. 

^radiator  =  ^recycle  —  ^productrecycled  (21) 


^productrecycled 


^feedQ^feed  (^target  ^ambient)  “E  ^recycle Q^product  (^target  ^radiator) 

Q^product  (^product  —  ^radiator^) 


(20) 


exchanger  tubes.  The  parameters  used  in  the  radiator  models  are 
shown  in  Table  3.  The  fin  is  constructed  from  aluminum  and  the 
tube  from  stainless  steel. 

As  with  the  reactor  model,  Simulink  solves  the  temporal  com¬ 
ponents  of  the  equations  while  the  spatial  components  are  solved 
by  creating  a  series  of  10  differential  elements. 

3.4.  System  controller 

The  amount  of  material  that  needs  to  be  recycled,  in  addition  to 
the  material  that  needs  to  be  pumped  to  the  reactor,  is  regulated  by 
a  control  system.  The  overarching  goal  of  the  controller  is  to 
regulate  the  flow  rate  of  fresh  slurry  into  the  reactor  based  on 
requested  hydrogen  from  the  fuel  cell  and  ballast  tank  pressure. 


In  this  way,  both  constraints  for  a  sufficient  preheat  and 
adequate  dilution  for  maximum  temperature  maintenance  could 
be  met.  The  feed  temperature  target  temperature  is  120  °C  and  the 
product  maximum  temperature  is  250  °C. 

With  the  endothermic  reaction  of  alane,  heat  needs  to  be  added 
to  the  reactor  continuously.  As  a  result,  a  slightly  different 
controller  for  heat  flux  to  the  reactor  exterior  surface  was  imple¬ 
mented.  Rather  than  just  being  a  constant  used  only  during  startup, 
as  is  the  case  for  the  AB  reactor,  the  heater  for  the  alane  system 
controller  is  based  on  the  hydrogen  requested  and  the  heat  of  re¬ 
action  of  the  storage  material  and  the  heat  required  to  get  to  the 
target  temperature  in  the  form  of  Eq.  (22). 


Table  3 

Parameters  used  in  the  heat  exchanger  models. 


Slurry  parameter 

Value 

Units 

Gas  parameter 

Value 

Units 

Fin  diameter  (dfin) 

0.025 

m 

Fin  diameter  (dfin) 

0.016 

m 

Heat  exchanger  tube  diameter  (dhxoo) 

0.0095 

m 

Heat  exchanger  tube  diameter  (dhxoo) 

0.00635 

m 

Fin  thickness 

0.00046 

m 

Fin  thickness 

0.00046 

m 

Fin  pitch  (pfin) 

0.00254 

m 

Fin  pitch  (pfin) 

0.00254 

m 

Twist  ratio  (y  =  diameters/1 80°  twist) 

9.84 

Twist  ratio  (y  =  diameters/1 80°  twist) 

9.84 

Slurry  HX  length  AB 

1.6 

m 

Gas  HX  length  AB 

1.25 

m 

Slurry  HX  length  alane 

1.2 

m 

Gas  HX  length  alane 

1.0 

m 
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eater 


^H2, required  (i^l  ^^rxn,l  T"  /^2^H  rxn,l ) 

(01  +  02 ) 

+  /TlfgedCp  feed  (^target  —  Tfeed)  +  Qpi  (^reactor  —  ^target) 


(22) 

To  maintain  the  reactor  temperature  near  the  target  value,  the 
heat  is  adjusted  up  or  down  by  a  PI  controller  based  on  the  differ¬ 
ence  between  the  reactor  and  target  temperature. 


3.5.  Recuperator  for  alane  slurry 

In  the  case  of  alane,  a  recuperator  is  used  to  provide  heat  to  the 
reactor.  The  recuperator  is  assumed  to  be  a  counterflow  tube-in¬ 
tube  heat  exchanger.  It  is  modeled  by  a  series  of  four  differential 
equations  representing  the  slurry  feed  (Tfeed)  and  product  tem¬ 
peratures  (Tproduct)  and  inner  (Tinner)  and  outer  (Touter)  wall  tem¬ 
peratures.  Heat  convection  is  accounted  for  as  if  flows  from  the  hot 
product  to  the  inner  wall  to  the  cooler  product  using  Eq.  (23). 


Table  5 

Parasitic  power  requirements. 


Component 

Value/equation  (Watts) 

Operation 

Volume  displacement 

20 

Start-up  only 

tank  mixers 

Pumps 

5.04  x  105  AP  (bar)  * 

Variable 

^slurry  /Pslurry 

Radiator  blower 

91 

Constant 

^H2,  total 


at 


nH2,in  -  nH2,FC 


(25) 


To  maintain  hydrogen  pressure  above  5  bar  during  startup,  the 
ballast  tank  volume  was  required  to  be  30  L. 


3.7.  Parasitic  power 

The  parasitic  power  includes  system  pumps,  heaters,  blowers, 
and  mixers  that  require  electricity  either  continuously  or  periodi¬ 
cally  depending  on  system  operations  (  Table  5).  This  parasitic  po- 


AnnerPfeedQ>  (  — +  ufeed 


mner^inner71'  (  ^feed  ^inner  )  —  0 


^Q^inner^ 


Annerwall Pmetal^-p, metal  (  ^inner^inner7r  (  ^inner 


^feed  )  3“  ^innerwall^innerwallTC: (  ^inner  ^product)  —  ^ 


(QJ  QJ  \ 

T  ^product  ^  J  "T  ^innerwall^innerwallTC:  (^product  —  ^inner)  3“  ^outerf^outer71' ^product  —  ^outer ^  —  0 


^outerPmetal ^p, metal  ^  ^  3"  frouterwall^outerwall7^ 


outer  ^product )  —  ^ 


(23) 


In  these  equations,  the  cross  sectional  areas  04)  and  diameters 
(D)  represent  the  inner  tube  (tube  center)  and  tube  annulus  as  well 
as  the  inner  and  outer  of  the  metal  tubes  ( Table  4).  The  heat  transfer 
coefficient  for  the  feed  is  based  on  the  laminar,  fully  developed  flow 
Nusselt  number  of  4.36  for  constant  heat  flux.  The  heat  transfer 
coefficients  for  the  annulus  are  provided  in  Eq.  (9)  and  Eq.  (10).  As 
with  the  reactor  and  the  radiators,  Simulink  solves  the  temporal 
components  of  the  equations  while  the  spatial  components  are 
solved  by  creating  a  series  of  10  differential  elements. 


wer  is  calculated  and  summed  throughout  the  model  run  and  is 
used  to  determine  the  storage  system  onboard  efficiency. 

The  model  also  predicts  the  storage  system  mass  and  volume 
based  on  user-inputted  dimensions  such  as  mass  of  the  CH,  radiator 
length,  ballast  tank  volume,  and  the  reactor  design.  These  values 
were  estimated  based  on  simple  system  designs  and  the  weight 
and  volume  of  commercial  components.  They  are  provided  in 
Table  6.  The  total  system  mass  and  volume  can  then  be  used  during 
tradeoff  studies  to  study  impacts  of  the  system  parameters. 


3.6.  Ballast  tank 


The  pressure  of  the  ballast  tank  varies  based  on  the  tank  tem¬ 
perature  and  the  number  of  moles  in  the  tank  assuming  ideal  gas 
law.  The  temperature  is  calculated  based  on  the  temperature  of  the 
hot  gases  entering  the  ballast  tank  from  the  radiator  using  Eq.  (24). 

^H2, total Cp,H2  ^  ^  +  ^H2, in Qj,H2  (^radiator  —  ^Bt)  —  0  (24) 

The  number  of  moles  in  the  tank  is  based  on  the  number  of 
moles  entering  the  tank  from  the  reactor  and  the  number  of  moles 
leaving  to  meet  the  drive  cycle  as  shown  in  Eq.  (25). 


Table  4 

Parameters  used  in  the  recuperator  model. 


Slurry  parameter 

Value 

Units 

Inner  tube  diameter  (Dinnerwaii) 

0.0127 

m 

Outer  tube  diameter  (Douterwaii) 

0.0196 

m 

Wall  thickness 

0.00056 

m 

Table  6 

Parameters  used  to  calculate  the  overall  system  mass  and  volume. 


Component 

Mass 

coefficient 

(kg) 

Volume 

coefficient 

(L) 

Determining  dimension 

CH  material 

M 

0 

M  =  mass  of  CH,  kg 

Volume 

exchange 

51.9  *  M/p  +  1.75 

M/p  *  1000 

M  =  mass  of  CH,  kg 
p  =  slurry  density,  kg  m  3 

tank 

Reactor 

Based  on  stainless  steel  and  dimensions 

Gas  radiator 

0.32  *  L 

0.2  *  L 

L  =  radiator  length,  m 

Slurry  radiator 

0.54  *  L 

0.87  *  L 

L  =  radiator  length,  m 

Ballast  tank 

0.17  *  V 

V 

V  =  tank  volume,  liters 

Ammonia 

0.0032  *  C 

0.0068  *  C 

C  =  concentration,  ppm 

scrubber 

Borazine 

0.0017  *  C 

0.0015  *  C 

C  =  concentration,  ppm 

scrubber 

Recuperator 

Based  on  stainless  steel  and  dimensions 

Balance  of  plant 

32.1 

26.8 
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4.  Model  validation  experiments 

Key  components  were  tested  in  the  laboratory  individually  and 
their  performance  compared  to  their  respective  models  as  a  means 
of  validating  the  overall  system  model. 


4.2.  Reactor  testing 

Once  the  model  was  developed,  the  reactor  piece  of  the  model 
developed  as  described  above  was  then  compared  to  the  pre¬ 
liminary  experimental  results  obtained  by  Los  Alamos  National 
Laboratory  (LANL)  using  a  prototypical  flow  through  reactor  scaled 
down  from  the  design  anticipated  for  automotive  application.  The 
reactor  design  consisted  of  a  1-cm  inner  diameter  tube  with  a  19.1- 
cm  long  heated  section.  A  McMaster  Carr  deep  hole  drip  bit  rotated 
at  40  rpm  inside  of  the  tube  as  a  means  of  mixing  and  moving  the 
slurry  through  the  reactor  as  shown  in  Fig.  3.  The  reactor  temper¬ 
ature  was  controlled  on  the  outer  surface  over  a  range  between  150 
and  235  °C.  Alane  slurries  between  20  and  60  wt%  and  AB  slurries  at 
20  wt%  flowed  through  the  reactor.  The  auger  speed  was  adjusted 
between  12  and  40  rpm  and  the  residence  time  between  2.4  and 
7.6  min.  For  each  condition,  the  reactor  was  allowed  to  equilibrate 
and  the  hydrogen  generation  rate  was  measured  using  a  calibrated 
flow  meter.  The  hydrogen  flow  rate  then  was  compared  to  the  flows 
predicted  by  the  model. 

The  automotive  dimensions  and  flows  used  in  the  model  were 
scaled  down  to  those  used  in  the  LANL  experiment.  The  tempera¬ 
ture  within  the  small-scale  reactor  was  assumed  to  be  constant  at 
the  controlled  value.  Reaction  enthalpies  and  heat  losses  were  not 
included  as  assumed  in  automotive  case.  The  tanks,  heat  ex¬ 
changers,  and  other  balance-of-plant  components  used  in  the 
automotive  storage  system  model  were  removed  to  predict  the 
reactor  behavior.  Only  the  minimum  slurry  density  (psiurry,min)  was 
adjusted  to  fit  the  experimental  data.  The  value  of  34  kg  m-3  for 
alane  was  fit  to  the  data  in  a  least  squares  regression  analysis.  Using 
this  minimum  density  value,  the  coefficient  of  determination  was 
calculated  to  be  90%  for  12  data  points  fit  in  Table  7. 

The  data  include  a  range  of  solids  loadings,  residence  times, 
auger  speeds,  and  reaction  temperatures.  There  is  reasonably  good 
fit  between  the  model  and  the  experiment  over  the  range  of  con¬ 
ditions.  Using  the  parameters  obtained  from  the  data  set  modeled 
in  able  6,  a  separate  set  of  alane  data  were  fit  with  the  model. 
These  data  were  obtained  for  a  20  wt%  slurry,  residence  time  of 
6.8  min,  and  an  auger  speed  of  40  rpm  over  a  temperature  range  of 
150-275  °C.  The  results  of  this  data  fit  are  shown  in  Fig.  4.  While 
the  fit  is  reasonably  good,  the  model  does  tend  to  slightly  under¬ 
predict  hydrogen  generation  rates  at  the  low  temperatures  and 
over-predict  them  at  the  higher  temperatures. 

This  same  model  was  also  used  to  fit  the  experimental  data 
developed  for  AB.  Only  three  data  points  at  20  wt%  solids  were 


Feed 


Fig.  3.  Flow-through  reactor  design  used  to  validate  the  model. 


Table  7 

Comparison  of  measured  and  model  predicted  alane  conversion. 


Solids 

loading 

(wt%) 

Reactor 

residence 
time  (min) 

Auger 

speed 

(rpm) 

Average 

reaction 

temperature 

(°C) 

Measured 

alane 

conversion 
(mol  mol  :) 

Modeled 

conversion 
values 
(mol  mol-1) 

50% 

7.6 

12 

185 

16% 

11% 

50% 

7.6 

40 

185 

11% 

4.6% 

50% 

4.2 

12 

187 

7% 

6.6% 

50% 

7.6 

12 

214 

88% 

80% 

50% 

7.6 

40 

214 

74% 

53% 

50% 

4.2 

12 

214 

38% 

49% 

20% 

6.8 

40 

188 

10% 

7.1% 

20% 

6.8 

40 

212 

38% 

50% 

20% 

6.8 

40 

235 

84% 

100% 

60% 

7.2 

12 

180 

5% 

6.5% 

60% 

7.2 

12 

194 

20% 

21% 

60% 

7.2 

12 

208 

48% 

55% 

Average  Reaction  Temperature  (°C) 

Fig.  4.  Comparison  of  model  predicted  hydrogen  flow  rate  to  that  measured  experi¬ 
mentally  for  20  wt%  alane  Slurry. 

measured  due  to  difficulties  in  pumping  the  AB  slurry.  In  this  case, 
the  density  of  the  slurry  product  of  370  kg  m  3  was  measured  and 
used  in  the  model  as  the  minimum  slurry  density  (psiurry,min)-  The 
model  predicted  the  experiment  accurately  at  150  °C,  the 
maximum  temperature  of  the  batch  data,  but  did  not  extrapolate 
well  to  higher  temperature  values  (Fig.  5).  As  a  result,  the  kinetic 


1000 

900  Model  Prediction 


Average  Reaction  Temperature  (°C) 


Fig.  5.  Comparison  of  model  predicted  hydrogen  flow  rate  to  that  measured  experi¬ 
mentally  for  20  wt%  AB  Slurry. 


512 


K.P.  Brooks  et  al.  /  Journal  of  Power  Sources  271  (2014)  504—515 


Table  8 

Updated  kinetic  parameters  based  on  measured  AB  values. 


Parameter 

Units 

AB 

1st  equivalent 

2nd/3rd  equivalent 

Arrhenius  factor  (Ai) 
Activation  energy  (Ei) 

s_1 

kj  mol-1 

2.81  x  102 

29.9 

4.18  x  1015 

149 

parameters  were  adjusted  to  fit  these  three  experimental  data 
points.  These  updated  kinetic  parameters  are  shown  in  Table  8  and 
are  used  for  all  subsequent  analysis. 


4.2.  Radiator  testing 

Validation  of  the  radiator  models  was  performed  using  eight 
finned  heat  exchangers  in  series  assembled  into  a  rectangular 
configuration  as  shown  in  Fig.  6.  These  finned  tubes  produced  by 
Energy  Transfer  MDE  were  made  from  0.95  cm  diameter  stainless 
steel  tubes  with  a  total  length  of  2.43  m.  The  aluminum  fins  brazed 
onto  the  tube  were  0.76  cm  tall  and  0.046  cm  thick,  and  were 
spaced  5.5  fins  cm-1.  Inside  the  tube  was  a  stainless  steel  turbulator 
with  a  pitch  of  7.6  cm  per  revolution.  Ambient  air  at  ~29  °C  was 
blown  through  this  radiator  at  a  constant  rate  of  6  m  s_1. 

Because  of  the  cost  of  both  the  raw  materials  and  the  prepara¬ 
tion  of  CH  slurries,  a  simulant  was  developed  to  be  used  in  the 
radiator  test  to  match  that  of  an  AB  slurry  instead  of  using  reacted 
CH  materials.  A  40  wt%  polyimide  (P84-N,  10  to  59  microns)  slurry 
in  the  AR-20  silicone  oil  was  selected  as  a  good  physical  simulant, 
and  it  modeled  a  spent  40  wt%  AB  slurry  reasonably  well.  The 
simulant  slurry  concentration  was  selected  because  it  matched  the 
viscosity  of  the  product— both  were  -200  cP  with  a  shear  stress  of 
1  Pa.  The  simulant  slurry  had  a  higher  wt%  and  vol%  solids  than  the 
AB  slurry,  but  it  bound  both  of  these  values.  The  two-phase  system 
was  selected  to  tolerate  the  elevated  temperatures  and  pressures 
required  for  the  validation  testing. 

The  polyimide  slurry  was  electrically  heated  in  a  stirred  feed 
vessel  to  between  150  and  230  °C.  Once  the  set  point  was  reached, 
the  slurry  was  pumped  through  the  radiator  at  between  44  and 
263  g  min-1,  and  the  inlet  and  outlet  temperatures  across  the  heat 
exchanger  were  measured.  The  model  results  were  then  compared 
to  the  experimental  results  as  shown  in  Fig.  7.  The  good  fit  between 
the  model  and  experimental  indicates  that  the  model  is  adequately 
validated. 


Fig.  6.  Prototypical  heat  exchanger  used  to  validate  the  model. 


Fig.  7.  Comparison  of  the  experimental  and  modeling  results  for  the  heat  exchanger. 

4.3.  Validation  of  other  components 

In  addition  to  the  reactor  and  radiator,  other  components  were 
validated  as  part  of  this  project.  These  components  include  the 
volume  displacement  tank,  the  phase  separator,  and  the  impurities 
clean-up  system.  These  components  require  validation  as  part  of 
the  system  development,  but  complete  phase  separation  and  im¬ 
purities  clean-up  is  assumed  in  the  model.  As  a  result  their  vali¬ 
dation  is  not  critical  to  the  modeling  effort  and  is  provided 
elsewhere  [18  . 

5.  Simulation  results 

A  variable  step  solver  “ode23tb”  is  used  in  Simulink  to  solve  the 
model  for  all  studies  performed  in  this  work.  To  avoid  algebraic 
errors  in  the  model  resulting  from  reactor  outputs  feeding  back  to 
the  controller,  unit  delays  are  used.  To  demonstrate  the  usefulness 
of  the  simulation  for  real  vehicle  driving  scenarios,  transient  sim¬ 
ulations  using  various  drive  cycles  were  run.  The  test  matrix  was 
provided  by  our  HSECoE  partners  to  test  various  storage  system 
models  developed  by  team  partners  [17]  and  includes  1)  a  highway 
cycle  (UDDS  +  HWFET),  2)  an  aggressive  driving  cycle  (US06),  3)  a 
cold-start  cycle  (cold  FTP),  and  4)  a  hot  cycle  (SC03-Hot).  The 
hydrogen  demand  data  from  the  fuel  cell  was  provided  by  the 


Time  (sec)  x  104 

Fig.  8.  Ballast  tank  pressure  as  a  function  of  time  for  50  wt%  AB  with  the  US06  Drive 
Cycle. 


K.P.  Brooks  et  al.  /  Journal  of  Power  Sources  271  (2014)  504—515 


513 


National  Renewable  Energy  Laboratory,  an  HSECoE  partner  [23  .  All 
drive-cycle  simulations  were  successfully  implemented  into  the 
model  but  for  conciseness,  only  important  simulation  results  are 
provided  in  this  paper.  All  the  drive  cycles  were  repeated  until  the 
DOE-required  5.6  kg  of  usable  hydrogen  was  produced. 

The  US06  drive  cycle,  the  most  aggressive  drive  cycle  considered 
in  this  study,  was  used  as  the  hydrogen  requirement  input  to  the 
model.  As  is  consistent  with  all  drive  cycles,  the  pressure  in  the 
ballast  tank  shows  an  initial  decrease  with  no  hydrogen  being 
generated  while  the  reactor  is  being  heated  to  temperature  as 
shown  in  Fig.  8.  This  first  ~2-3  min  of  the  US06  drive  cycle  de¬ 
termines  the  size  of  the  ballast  tank.  The  tank  must  be  large  enough 
to  maintain  the  hydrogen  pressure  above  5  atm  during  startup. 
Once  the  reaction  is  initiated,  the  hydrogen  pressure  in  the  ballast 
tank  varies  only  slightly  around  the  25  atm  set  point  until  the  slurry 
is  consumed  and  the  ballast  tank  pressure  is  lost. 

Although  the  US06  drive  cycle  reaches  hydrogen  usage  of  up  to 
0.8  mol  H2/second  (80  kWe),  the  storage  system  ballast  tank  pro¬ 
vides  a  buffer.  As  a  result,  with  a  50  wt%  AB  slurry  for  example,  the 
storage  system  is  required  to  produce  a  maximum  of  only  0.43  mol 
H2/second  (43  kWe)  as  shown  in  Fig.  9.  The  alane  system  is  required 
to  produce  a  maximum  of  0.32  mol  H2/second  (32  kWe).  This 
significantly  reduces  the  system's  size  and  capacity  requirements. 
During  startup,  an  electric  heater  brings  the  reactor  to  temperature 
before  the  reaction  is  initiated.  During  the  entire  operation,  electric 
pumps  and  mixers  are  used  to  maintain  the  flow  and  ensure  ho¬ 
mogeneity  of  the  slurry  mixture.  All  of  these  electrical  components 
result  in  parasitic  loss  of  hydrogen  that  cannot  be  used  for  primary 
automotive  propulsion.  The  DOE  has  given  the  requirement  for  90% 
onboard  efficiency.  Based  on  the  amount  of  electrical  power 
required  and  using  0.02  g  H2/s/kWe  conversion  of  hydrogen  to 
electrical  power,  the  onboard  efficiency  was  calculated  using  Eq. 
(26): 


Efficiency 


^H2, total  ^H2,parasitics 

^H2,  total 


(26) 


where  mH2  total  is  the  mass  of  hydrogen  generated  by  the  system 
and  mH2  parasitics  is  the  amount  of  hydrogen  that  must  be  converted 
to  electrical  power  for  parasitic  losses.  For  the  US06  drive  cycle  as 
well  all  other  drive  cycles  studied  with  50  wt%  AB,  the  onboard 


efficiency  is  predicted  to  be  >95%.  In  contrast,  the  onboard  effi¬ 
ciencies  for  50%  alane  are  predicted  to  be  87-93%.  The  alane  re¬ 
action  is  endothermic,  and  as  a  result,  it  requires  reactor  heating 
during  the  entire  operations.  In  contrast,  AB  produces  enough  heat 
that  no  reactor  heating  is  required.  Reactor  heating  is  the  major 
parasitic  loss  in  the  alane  system.  Such  components  as  the  volume 
displacement  tank,  phase  separator,  and  impurities  clean-up  sys¬ 
tem  do  not  impact  the  on-board  efficiency  of  the  system  since  they 
do  not  require  power  for  their  operation. 

The  high  hydrogen  requirement  of  the  US06  also  results  in  the 
most  stringent  case  for  the  storage  system  radiators.  As  shown  in 
Fig.  1,  there  is  one  radiator  to  cool  the  separated  hydrogen  and  one 
for  the  spent  slurry.  The  DOE  targets  require  that  the  hydrogen 
going  to  the  fuel  cell  be  maintained  at  a  temperature  <85  °C.  There 
is  not  a  specific  target  for  the  cooling  required  for  the  slurry,  but  at 
sustained  temperatures  >60  °C,  the  fresh  CH  that  is  in  contact  with 
the  spent  CH  through  only  a  membrane  in  the  volume  displace¬ 
ment  tank  will  begin  to  generate  hydrogen.  As  a  result,  the  slurry 
and  hydrogen  radiators  have  been  sized  to  maintain  their  exit 
temperature  less  than  the  required  values.  The  temperatures  of  the 
hydrogen  and  slurry  for  both  AB  and  alane  are  shown  in  Figs.  10  and 
11,  respectively.  As  can  be  seen  from  these  results,  startup  produces 
the  highest  hydrogen  exit  temperature  for  both  AB  and  alane  and 
for  the  AB  exit  slurry  temperature.  In  contrast  the  slurry  temper¬ 
ature  for  alane  increases  over  time. 

The  model  demonstrated  that  the  storage  system  as  designed 
could  meet  the  drive-cycle  requirements  for  both  alane  and  AB  for 
the  other  three  drive  cycles  as  well.  Similar  results  were  obtained 
for  the  other  drive  cycles  and  will  not  be  provided  here. 

The  system  mass  and  volume  also  were  calculated  as  part  of  this 
modeling  effort.  Based  on  the  current  designs,  the  AB  system  has  a 
mass  and  volume  of  126  kg  and  133  L,  respectively.  For  hydrogen 
production  of  5.6  kg,  the  gravimetric  and  volumetric  capacities 
calculated  for  the  AB  system  are  0.044  kg  H2/kg  system  and 
0.042  kg  H2/L  system.  The  alane  system,  in  contrast,  has  a  mass  and 
volume  of  172  kg  and  153  L,  respectively.  For  hydrogen  production 
of  5.6  kg,  the  gravimetric  and  volumetric  capacities  calculated  for 
the  alane  system  are  0.033  kg  H2/kg  system  and  0.037  kg  H2/L 
system.  The  primary  difference  in  the  mass  and  volume  between 
these  two  systems  is  the  additional  mass  and  volume  of  CH 
required  for  the  alane  system.  Neither  AB  nor  alane  systems  meet 
the  2017  DOE  technical  targets  for  gravimetric  and  volumetric 
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Fig.  9.  Comparison  of  hydrogen  required  by  the  drive  cycle  and  hydrogen  produced  by  Fig.  10.  Temperature  of  gas  and  slurry  radiator  products  for  50  wt%  AB  with  the  US06 
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Fig.  11.  Temperature  of  gas  and  slurry  radiator  products  for  50  wt%  alane  with  the 
US06  Drive  Cycle. 
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Fig.  13.  Sensitivity  analysis  comparing  the  slurry  storage  capacity  to  the  predicted 
system  gravimetric  capacity  for  both  an  endothermic  and  an  exothermic  system. 


capacities  of  0.055  kg  Fh/kg  system  and  0.04  kg  H2/L  system. 
Additional  evaluation  of  these  systems  relative  to  the  technical 
targets  is  provided  elsewhere  [18  . 

6.  Sensitivity  analysis 

The  DOE  does  not  assume  that  AB  or  alane  slurries  will  be  the 
materials  ultimately  used  in  vehicles.  It  is  assumed  that  additional 
chemical  hydrogen  storage  material  discovery  and  development 
will  be  performed.  As  a  result,  one  of  the  purposes  of  the  model  is  to 
help  guide  potential  material  developers  regarding  the  range  of 
material  properties  that  could  be  implemented  successfully  into  a 
chemical  storage  system  similar  to  what  was  developed  here.  To 
this  end,  the  model  also  was  used  to  perform  a  sensitivity  analysis 
and  better  understand  the  impact  of  the  material  properties  being 
able  to  meet  the  DOE  targets.  The  models  were  run  for  3000  s  using 
five  US06  drive  cycles  back-to-back.  Key  material  properties  were 
varied  individually  around  the  values  used  in  the  model  described 
above  until  one  of  the  DOE  targets  was  not  met  or  the  model  could 
not  run  successfully. 

Key  parameters  of  a  hydrogen  storage  material  are  heat  of  re¬ 
action,  hydrogen  capacity,  and  reaction  rate.  These  parameters 
were  explored  for  both  alane  and  AB  as  surrogates  for  possible 


Fig.  12.  Sensitivity  analysis  comparing  heat  of  reaction  to  the  predicted  onboard  effi¬ 
ciency  for  both  an  endothermic  and  an  exothermic  system. 


future  hydrogen  storage  materials.  The  heat  of  reaction  was  varied 
for  these  two  materials  and  the  results  were  evaluated.  For  either 
exothermic  or  endothermic  materials,  as  the  heat  of  reaction  moves 
further  away  from  zero,  the  reactor  temperature  increases.  This 
results  in  higher  ballast  tank  and  product  tank  temperatures.  At 
some  point,  temperatures  may  become  excessive.  In  addition,  as 
the  reaction  heat  becomes  increasingly  exothermic,  the  onboard 
efficiency  and  conversion  decrease  as  well.  Based  on  the  results  of 
the  model,  onboard  efficiency  drops  below  the  DOE  technical  target 
value  of  90%  at  -11  kj  mol-1  H2  as  shown  in  Fig.  12. 

The  hydrogen  capacity  of  the  material  impacts  the  onboard  ef¬ 
ficiency  and  gravimetric  and  volumetric  capacity.  Lower  hydrogen 
capacity  decreases  the  onboard  efficiency,  similar  to  the  heat  of 
reaction.  Values  of  -3.75%  hydrogen  for  an  alane  slurry  and  2.7  wt% 
hydrogen  for  an  AB  slurry  result  in  <90%  onboard  efficiency.  More 
significantly,  the  gravimetric  density  is  decreased  with  low 
hydrogen  capacity.  For  both  AB  and  alane,  the  DOE  technical  target 
gravimetric  factor  of  0.055  can  be  met  at  hydrogen  capacities  of 
>11.4%  (See  Fig.  13).  The  system  mass  was  similar  for  both  systems 
in  spite  of  their  differences.  For  a  50  wt%  slurry,  this  would  require 
125%  and  50%  increases  in  hydrogen  capacity  for  alane  and  AB, 
respectively,  relative  to  their  current  values.  Previous  sensitivity 
analyses  indicate  that  the  volumetric  capacity  will  be  met  when  the 
gravimetric  capacity  is  achieved. 

An  increase  in  reaction  kinetics  above  their  current  values  has 
little  effect  on  the  drive-cycle  performance  or  the  system  param¬ 
eters.  Similarly,  decreasing  the  reaction  kinetics  has  little  impact  as 
well  up  to  the  point  that  conversion  becomes  so  poor  that  the  drive 
cycle  cannot  be  met.  This  point  occurs  at  -50%  of  the  reaction  rate  of 
AB  and  -2%  of  the  reaction  rate  of  alane. 

7.  Conclusions 

To  determine  the  feasibility  of  chemical  hydride  materials  for 
onboard  hydrogen  storage  for  transportation  based  fuel  cell  appli¬ 
cations,  representative  models  were  developed  for  exothermic  and 
endothermic  slurry  systems.  Ammonia  borane  and  alane  were 
chosen  as  the  materials  for  proof-of-concept  studies,  and  the  re¬ 
sults  were  compared  to  experimental  data  for  a  flow-through 
reactor  system  and  a  prototypical  radiator.  The  fits  between  the 
model  results  and  the  experimental  results  were  good  for  both 
components,  providing  validation  of  the  models. 
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Once  validated,  these  models  were  evaluated  relative  to  the  DOE 
technical  targets  and  a  set  of  four  drive  cycles.  Preliminary  obser¬ 
vations  from  the  simulations  include  1 )  the  drive  cycles  could  be 
met  for  both  the  exothermic  and  endothermic  models,  2)  both 
systems  will  require  higher  hydrogen  capacity  than  is  possible  for  a 
50  wt%  slurry  to  meet  the  DOE  gravimetric  capacity  target  of 
0.055  kg  Eh/kg  system,  and  3)  the  DOE  onboard  efficiency  target  of 
90%  can  be  met  for  all  exothermic  reaction  enthalpies  and  for 
endothermic  reaction  enthalpies  <11  kj  mol-1  H2. 
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Nomenclature 

A:  area  (m2) 

A':  Arrhenius  factor  (s-1) 

C;  concentration  (ppm) 

Cp ;  heat  capacity  (J  kg-1  K) 


D:  radiator  diameter  (m) 
d:  reactor  diameter  (m) 

Ea:  activation  energy  (J  mol-1) 
f:  mass  fraction 

h:  heat  transfer  coefficient  (W  m-2  K) 
k:  reaction  rate  (s-1) 

/<';  thermal  conductivity  (W  m-1  K) 
k':  constant  for  (J  calculation  (1/°C) 

L:  length  (m) 

m;  mass  flow  rate  (kg  s-1) 

M:  mass  of  CH  (kg) 

MW:  molecular  weight  of  chemical  hydride 

n;  molar  flow  rate  (mol  s-1) 

n:  exponent  /moles  of  hydrogen 

Nu:  Nusselt  number 

P;  pressure  (atm) 

pfin:  fin  pitch 

Pr:  Prandtl  Number 

Q;  heat  (kW  m-1) 

r*:  ratio  of  inner  and  outer  diameter  of  slurry 
R:  universal  gas  constant  (J  mol-1  K) 

Re:  Reynolds  Number 
T:  temperature  (K) 
t:  time  (s) 
u:  velocity  (ms-1) 

V:  volume  (L) 
x:  mass  fraction 
y:  turbulator  twist  ratio 
z:  axial  distance  (m) 

Z;  axial  displacement  per  revolution  (m) 
a:  extent  of  reaction 

/?:  stoichiometric  ratio  of  hydrogen  to  CH 
/?';  constant  for  calculating  (3  as  a  function  of  T 
5fin:  surface  area  ratio  for  the  radiator  fins 
AH:  reaction  enthalpy  (kj  mol-1  H2) 
p:  viscosity  (N-s  m-2) 

77;  fin  efficiency 
v:  rotational  speed  (rpm) 
p:  density  (kg  m-3) 

Subscripts 

air:  relating  to  air 

ambient:  ambient  temperature 

ballast:  associated  with  the  ballast  tank 

CH:  chemical  hydride 

CHP:  chemical  hydride  product 

exterior:  outer  surface  of  the  radiator  tube 

Hzfc -  hydrogen  to  fuel  cell 

H^in-  hydrogen  into  device 

H2,  paras  tic-  amount  of  hydrogen  required  to  power  parasitic  equipment 

H2, total-  total  hydrogen  available 

feed:  feed 

fin:  fin  dimensions 

Hy.  hydrogen  gas 

heater:  heater 

i:  reaction  number  (e.g.  1,2) 

in:  inner  diameter  reactor 

inner:  inner  diameter  of  recuperator  inner  tube 

inert:  inert  material  (slurrying  agent) 

innerwall:  outer  wall  of  recuperator  inner  tube 

interior:  inner  surface  of  the  radiator  tube 

loss:  heat  loss 

metal:  metal  for  reactor  or  radiator 

min:  minimum  allowable  value 

out:  outer  diameter  reactor 

outer:  inner  wall  of  the  outer  recuperator  tube 

outerwall:  outer  wall  of  the  outer  recuperator  tube 

PP.  proportional-integral  controller  value 

product:  associated  with  product 

product  recycle:  associated  with  cold  product  from  radiator  recycled  back  to  reactor 

radiator:  associated  with  the  radiator 

reactor:  associated  with  the  reactor 

required:  hydrogen  requested  from  the  fuel  cell 

recycle:  associated  with  hot  product  recycle  back  to  reactor 

rxn:  reaction  property 

set  point:  associated  with  a  set  point  value 

slurry:  slurry 

stir:  reactor  stirrer 

target:  associated  with  a  target  value 


